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Discontinuous quantum phase transitions besides their general interest are clearly relevant to 
the study of heavy fermions and magnetic transition metal compounds. Recent results show that 
in many systems belonging to these classes of materials, the magnetic transition changes from 
second order to first order as they approach the quantum critical point (QCP). We investigate here 
some mechanisms that may be responsible for this change. Specifically the coupling of the order 
parameter to soft modes and the competition between different types of order near the QCP. For 
weak first order quantum phase transitions general results are obtained. In particular we describe 
the thermodynamic behavior at this transition when it is approached from finite temperatures. This 
is the discontinuous equivalent of the non-Fermi liquid trajectory close to a conventional QCP in a 
heavy fermion material. 
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I. INTRODUCTION 

Recently, the subject of discontinuous quantum phase 
transitions has received much attention 1 . Besides 
the theoretical interest in these transitions there are 
many new experimental results showing that they oc- 
cur in heavy fermions and magnetic transition metal 
compounds 2,3 . In these materials as one approaches a 
quantum critical point, the transitions change their na- 
ture from second to first order—. There are many mech- 
anisms which can drive a second order transition into 
a first order one. For example, in compressible magnets, 
pressure may produce this changed. In antiferromagnets 4 
and superconductors^ this may be accomplished by an 
external uniform magnetic field . Here we will be inter- 
ested in more subtle mechanisms which arise due to a 
coupling of the order parameter to fluctuations^. In this 
case the effects are weaker being associated with what 
is generally known as weak first order transitions. These 
mechanisms become more relevant at very low tempera- 
tures where quantum fluctuations are important. 

The fluctuations we consider are of different types. In 
the first case they are soft modes and as a paradigm of 
this case we study the coupling of the superconductor or- 
der parameter to the electromagnetic field. This has been 
investigated by Halperin, Lubensky and Ma— but, since 
we are interested in the fully quantum mechanical case at 
zero temperature, our problem becomes similar to that 
treated by Coleman and Weinberg in particle physics^. 
Next, we consider the effect of coupling the order param- 
eter of a given phase to fluctuations of a different phase 
which competes with the former in the same region of 
the phase diagram. This problem is particularly rele- 
vant for heavy fermion systems where an inhomogeneous 
state consisting of antiferromagnetic and superconduct- 
ing regions has been observed^. This coupling between 
competing order parameters gives rise to new effects com- 
pared with the previous situation. In this case not only 
the nature of the quantum transitions may change but 
the ordered phase itself may become an inhomogeneous 
mixture with different values for the order parameter. 



We obtain general features of the behavior of a system 
approaching a weak first order quantum phase transition 
(WFOQPT) from finite temperatures. In spite that the 
correlation length does not diverge as T — > 0, scaling 
does apply for these WFOQPT and allows for a general 
description of this phenomenon. 



II. COUPLING TO SOFT MODES: THE 
COLEMAN- WEINBERG MECHANISM 

A well known case of quantum fluctuations inducing 
a weak first order transition is the Coleman -Weinberg 
mechanism 8 . In the solid-state version, we consider a 
superconductor represented by a complex field (^1,(^2) 
coupled to the electromagnetic fieldi2*ii. The Lagrangian 
density of the model is given by, 
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We are using H = c = 1 units and the indices fj,, v run 
from to d = 3. In Eq. (JJJ space and time are isotropic 
and consequently the dynamic critical exponent z = 1. 
For a neutral superfluid (q = 0) the system decouples 
from the electromagnetic field and has a continuous, zero 
temperature superfluid-insulator transition at m 2 = 
(see Fig. Pi . 

The effective potential method 1213-14 yields the quan- 
tum corrections to the action given by the Lagrangian 
density of Eq. (JJJ. At T = in the one loop approx- 
imation, the effective potential close to the transition 
(to w 0) is given by^> 
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where M comes from renormalization and is completely 
arbitrary^. We can take the value of M as the minimum 
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of the effective potential (ip). In this case we have the 
equation, for M = (ip), 
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which yields 



or alternatively for A 

A = 
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We can use the equation above to remove A from the 
effective potential Eq. (01 introducing another free pa- 
rameter, namely the non-trivial minimum of the effective 
potential {ip). In this case, Eq. reads 
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The change of the parameters above known as dimen- 
sional transmutation 8 can only be done if the effective 
potential has a non-trivial minimum {(p) 7^ as is clear 
from the equations above. Since close to the transition 
to s» 0, we conclude from Eqs. 10 and J3J that the con- 
dition for such a minimum to exist is 



A ~ q ■ 

Now let's define the coherence length 



and the London penetration depth 
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as is usually done for Ginzburg-Landau models. We can 
show that the condition of Eq. JJJ) is equivalent to 
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This is the same condition for the occurrence of weak first 
order transitions obtained in Ref. 7 . However, in general, 
if Xl ~ £, fluctuations of the order parameter can desta- 
bilize the Coleman- Weinberg mechanism and the result- 
ing quantum transitions can be second order (see Ref^). 

The condition (|1C)|> can be obtained from our results as 
follows. The ratio A_l/£, from Eq. © and Eq. @, is 



Ai = / \m 2 \ 
Substituting the result of Eq. (0J we get 
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This ratio is independent of to 2 as expected and it is 
clear that, if A ~ q 4 , we have (A L /£) < 1 or Ai C ^ 
as expected. We will discuss more about this condition 
later in section Ivl 

Let us consider that A ~ q 4 (or equivalently Xl <C £) 
and study the minima of the effective potential, Eq. JSJ. 
We find that at a critical value of the mass, given by 
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there is a first order transition to a superconducting 
state. Notice that the transition in the neutral super- 
fluid (q = 0) is continuous rather than first order and 
takes place at m 2 = 0. Therefore, quantum fluctuations 
due to the coupling to soft modes, the photons, lead to 
symmetry breaking in the normal region of the phase di- 
agram of the neutral superfluid. This occurs close to the 
quantum critical point of the chargeless system extend- 
ing the region where a symmetry broken phase appears 
in the phase diagram. The transition occurs at a finite 
value of the mass m 2 and is first order. The shift of the 
critical mass depends on a power of the coupling of the 
order parameter to the soft modes, in the present case 
the charge of the Cooper pairs. 

Introducing a parameter g = m 2 — m 2 which measures 
the distance to the transition, we can write the effective 
potential at T = as, 



Veff = ^(ip) 2 \m 2 -m 
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with the exponent value a = I reflecting the fact the 
transition is first order 11 . The associated latent heat is 



L h = 



calculated 11 . 



\(<p) 2 ■ Spinodals points at T — can also be 
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FIG. 1: Phase diagram of a charged superfluid coupled to 
photons. For completeness we show also the critical line of 
the neutral superfluid. Along the trajectory m 2 — mj? one 
can distinguish different regimes as explained in the text. 
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The finite temperature case can be studied replacing 
the frequency integrations in the calculation of the ef- 
fective potential by sums over Matsubara frequencies^. 
The generalization of the effective potential to finite tem- 
peratures close to the transition can be written as 11 



V eff (T) = - m 2 (yf\g\{l+ 
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and M(tp) = m + q f ■ The function Is(y) = I(y) can 
be obtained numerically integrating Eq. Ijl6|l . 

The finite temperature phase diagram is shown in 
Fig. ^ For completeness we show the critical line of the 
neutral superfluid, T$f = |w 2 |^, which is governed by 
the shift exponent ip -1 = z/(d + z — 2) = 1/2 in d = 3 
(see Refiifi). The new line of first order transitions is 
given by T c — yjm 2 



III. 



SCALING AT A WEAK FIRST ORDER 
QUANTUM TRANSITION 



We will now consider the system sitting at the new 
quantum phase transition point, i.e., at m 2 — to 2 and 
decrease the temperature. This is the equivalent of the 
non-Fermi liquid trajectory for ordinary quantum critical 
points in metallic systems. 




m/T 



FIG. 2: The temperature dependent specific heat along the 
trajectory m 2 = m 2 of the zero temperature first order tran- 
sition. The different regimes are discussed in the text. Notice 
the dashed straight line that indicates the InT behavior in 
regime II. 

For high temperatures, T 3> m c , which corresponds 
to the regime I of Fig. ^and Fig. |3 the function l3(y) 



saturates, I^(y < 0.12) ss —2.16. In this case the effective 
potential, 

1 f 4 32 T d+1 



and can be cast in the scaling form, 
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with F(0) — constant. This scaling form is reminiscent 
of that for the free energy close to a quantum critical 
point. In the present case of a discontinuous zero tem- 
perature transition, the critical exponent^ a — 1 and 
the characteristic temperature is, 
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with v = l/(d + z)ii. In this regime I or scaling regime, 



along the line to 



shown in Fig. 2] the free energy 



density has therefore the scaling form f(m — m c ,T) tx 
rp(d+z)/z anc L £ ne S p ec ifi c nea t; as shown in Fig.|2is given 

by, 



C/T\. 
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Then the thermodynamic behavior along the line m 2 — 
m 2 in regime I (T> m c ) is the same as when approach- 
ing the quantum critical point of the neutral superfluid, 
along the critical trajectory m 2 = 0. The system is un- 
aware of the change in the nature of the zero temperature 
transition and at such high temperatures charge is irrel- 
evant. 

When further decreasing temperature along the line 
m 2 = m 2 there is an intermediate, non-universal regime 
(regime II in Figs. ^ and |2J - In the present case for 
m c < T the specific heat C/T d ' z oc InT as can be seen 
from the dashed straight line shown in the semi-log plot 
of Fig. 

Finally, at very low temperatures, for T « m c and 
m 2 = to 2 , i.e., in regime III of Fig. the specific 
heat vanishes exponentially with temperature, C ' /T d / Z oc 
exp( — m c /T). The gap for thermal excitations is given by 
the shift m c of the quantum phase transition. The cor- 
relation length which grows along the line to 2 = to 2 with 
decreasing temperature reaches saturation in regime III 
at a value £s ~ m c 1 - Then, we can understand the expo- 
nential dependence of the specific heat as due to gapped 
excitations inside superconducting bubbles of finite size 
L ~ £. The gap between the states in the bubbles is 
A ~ L~ z ~ to c as we have found previously. 

Although the results above have been obtained for a 
particular model, the behavior in the scaling regime I and 
III should be universal and characteristic of any weak 
first order quantum transition. Notice that the relevant 
critical exponents which determine the scaling behavior 
in particular in regime I are those associated with the 
QCP of the uncoupled system which in the present case 
is the neutral superfluid. 
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IV. COUPLING BETWEEN ORDER 
PARAMETERS 

A. Quantum first order transitions in 
superconducting heavy-fermions 

Weak first order transitions and spontaneous symme- 
try breaking near quantum phase transitions can also oc- 
cur due to the coupling of an order parameter to nearly 
critical fluctuations of another phase. This provides a 
new mechanism for WFOQPT besides the coupling to 
gapless excitations as treated above and in Refi^ for ex- 
ample (see also Refsi 15 i 16 ). 

In this section we study the coupling, at T = 0, of 
an order parameter to fluctuations of a competing or- 
der. This type of coupling becomes important when two 
different phases are in competition on the same region 
of the phase diagram. In this case for some values of 
the microscopic interactions which are close, there are 
alternative ground states which interfere with one an- 
other. This is the situation, for example, of several heavy- 
fermionsii and of high-T c superconductors^ where in 
general the dispute is among superconductivity and an- 
tiferromagnetism. We obtain the effects on a given phase, 
of the coupling of its order parameter to fluctuations of 
the competing order, in the form of quantum corrections. 
We show that these can produce spontaneous symmetry 
breaking and change the nature of the transition even 
if these fluctuations are non critical. Again, due to the 
weak first order character of the new phase transition we 
expect to find scaling behavior as described previously. 

We consider the case of a heavy- fermion (HF) system 
with a superconducting phase that is close to an anti- 
ferromagnetic instability at T = as in Fig. [3J As 
before, we first describe the superconducting phase by 
the simplified Lorentz invariant free action as in Eq. 
but now coupled to an AF order parameter instead of the 
electromagnetic field. In many cases however the effects 
of dissipation should be taken into account. Then, later 
on we use another free action for the superconducting 
field. This different quadratic form is associated with a 
dynamic exponent z = 2, instead of z = 1 of the Lorentz 
invariant cas o 20 i 21 . 

Our model is of the Ginzburg-Landau type and con- 
tains three real fields. Two fields, 0i and 02, correspond 
to the two components of the superconductor order pa- 
rameter. The other field 03, for simplicity represents a 
one component antiferromagnetic order parameter. The 
results are immediately generalized to a three component 
AF vector field, with the unique consequence of changing 
some numerical factors^. The free functional associated 
with the magnetic part represented by the field 03, the 
sub-lattice magnetization, takes into account the dissipa- 
tive nature of the paramagnons near the magnetic phase 
transition^ in the metal and yields the propagator 




where r is a characteristic relaxation time and m 2 is re- 
lated to a local Coulomb repulsion U and the density of 
states at the Fermi level N(Ep) by 

m 2 p = 1 - UN(E F ). (19) 

The part of the action associated with the potential is 
given by 

y c/ (0i,0 2 ,0 3 ) = \rn 2 {4>\ + 4>l) + ^m 2 p 4>l+ 

2 ) + V p (<f> 3 ) + Vt{4n, 02, 3 ), (20) 
where the self-interaction of the superconductor field is 

K(0i,0 2 ) = ^(0 2 +0^) 2 , (21) 

and that of the antiferromagnet is given by 

Wa) = |0t (22) 

Finally, the last term is the (minimum) interaction be- 
tween the relevant fields, 

^(0i,0 2 ,0 3 )=u(0 2 + 0|)0i (23) 

This term is the first allowed by symmetry on a series 
expansion of the interaction. Notice that for u > 0, 
which is the case here, superconductivity and antiferro- 
magnetism are in competition and this term does not 
break any symmetry of the original model. However, in- 
cluding quantum fluctuations we show that spontaneous 
symmetry breaking can occur in the normal phase sepa- 
rating the SC and AF phases. This case is represented 
schematically in Fig. and is discussed in the subsequent 
sections. 




FIG. 3: Normal phase separating SC and AF phases. 



B. Superconducting transition 

We first consider the metal-superconductor transition 
in HF when the superconducting field is coupled, as in 
Eq. J23J), to AF fluctuations described by the propagator 
of Eq. (fTHjl fFig. Detailed calculation of the effective 
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potential have already been presented^. The general 
result is given by, 

V ef ^ c ) » iM^l + am^lM + ^t + O^ 5 ). (24) 

In Eq. (|24l) , A is a renormalized coupling but of the same 
order of the bare coupling A. The new coupling a intro- 
duced by fluctuations can produce a symmetry breaking 
in the normal state extending the SC region in the phase 
diagram at T — 0. The same mechanism turns the transi- 
tion to weak first order with a small latent heat^. Using 
the same argument as in the preceding section we can 
introduce again parameters associated with £ (coherence 
length) and Xl (London penetration depth) in this case 
given by 



(25) 
(26) 



It's difficult to carry out the dimensional transmutation 
analytically on the result of Eq. 1|24|) since the equations 
of minima are hard to solve. However, we can check nu- 
merically that the condition for the existence of minima 
away from the origin is equivalent to <C £ as in the 
previous case. It's also interesting to notice that the first 
order transition is produced by the cubic term in Eq. 124|) 
and this term is proportional to the magnetic mass m p . 
If magnetic fluctuations were critical, i.e., m p = 0, the 
only effect of the coupling would appear as a term pro- 
portional to (f> 5 . This term is usually neglected since its 
power is higher than those initially considered in the clas- 
sical potential and usually insufficient to create new min- 
ima around the origin. Therefore, if AF fluctuations were 
critical the effects of the quantum corrections in the tran- 
sition could be neglected. 



C. Antiferromagnetic transition 

In order to study the magnetic transition, as mentioned 
previously we consider two kinds of quadratic forms as- 
sociated with the free superconducting fields. The first is 
the usual Lorentz invariant free action that we have dis- 
cussed before in Sec. IIV Al Next we work with another 
free action which takes into account dissipation and is 
associated with a z — 2 dynamics. 

Let's first study the result for the Lorentz invariant 
propagator. Close to the AF transition we havei^ 



V. 



ef{93c) 



~2i4 



+ urn (j> 3c In 



(27) 



We notice directly from Eq. I|27|l that if the super- 
conductor order parameter fluctuations were critical we 
would obtain the same result of Eq. since in this case 



the parameter m 2 = 0. However, SC fluctuations are 
not critical but close to criticality and the last term of 
Eq. 127|) may become important. Of course, its relevance 
depends on the strength of the renormalized coupling u 
and the results considering this term leads to new and in- 
teresting changes in the ground state. We obtain, besides 
the two finite minima of the Coleman- Weinberg poten- 
tial of Eq. ©, two extra minima very close to the ori- 
giniS. The normal state associated with these minima 
has a small moment since the value of the order param- 
eter is related to the sub-lattice magnetization. An ad- 
ditional first order transition occurs when the other two 
minima away from the origin become the stable ones as 
the system moves away from the superconductor instabil- 
ity. This transition is from a small moment AF (SMAF) 
to a large moment AF (LMAF) and occurs even before 
the continuous mean field transition. When we move 
away from the magnetic transition the strength of this 
new term decreases with the value of m 2 and the two 
new minima move to the origin producing a normal state 
with vanishing magnetization again. It's important to 
note here that the SMAF is obtained because the mag- 
netic order parameter couples to fluctuations which are 
close to criticality but not critical. Critical fluctuations 
yield the same results of section ITT1 

Now, for many cases of interest, SC fluctuations are 
better described by a dissipative propagator associated 
with a z = 2 dynamica 20 i 21 similar to Eq. (JTSJ. This is 
useful to account for pair breaking interactions as mag- 
netic impurities that can destroy superconductivity 2 ^*^. 
It is given by, 



(28) 



The parameter m 2 is still related with the classical dis- 
tance from the phase transition and we have a relaxation 
time t'. In general we also have a non-dissipative term 
(as in the case studied above) but this term is neglected 
since the linear term is the most important in the low 
frequency region. Calculation of the effective potential is 
very similar to the previous cases and the result has the 
form 

Veff = \m 2 p 4>1 + + j^i^i + ™ 2 ) 5/2 ( 29 ) 

where the renormalized mass is 

2 



Mp = m ? ~ ~ 12^ 



m 3 p g 



and the renormalized coupling 

12 2 
g = g j mu ~ 



3 
8^ 



(30) 



(31) 



Quantum corrections can once again produce a weak first 
order transition. We note that in this case the renormal- 
ization of the potential obtained is easier since there are 
no logarithmic terms to renormalizeiS. Anyway, the re- 
sult is very similar to that presented in section llV Bl sincc 
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we have equivalent propagators for the fluctuations. In 
general the form of the quantum correction depends not 
only on the effective dimension of the uncoupled QCP 
but also on the form and specially the dynamics of the 
soft modes or competing fluctuations. 

The transformation (4>' 3 ) 2 = 2it</>| + m 2 makes the po- 
tential of Eq. i|29|) simpler and the analysis of its extrema 
shows that the transition can be first order depending on 
the coupling values. The appearance of SMAF phases 
is not possible in this case. We note, however, that the 
couplings are dependent on the renormalization param- 
eter M defined in Sec. ^ and a renormalization group 
approach is essential to make the results more reliable. 
In the next section we illustrate this with an analysis of 
the Coleman- Weinberg case. 



In this sense the effective potential is always given by 
the same equation with a suitable parametrization of the 
coupling. 

This argument can be formally stated in an renormal- 
ization group approach and we refer to the correspondent 
section of RefA As a result, a weak first order transi- 
tion occurs for any small values of A and q, since we 
can appropriately choose the renormalization mass. The 
only restriction comes from perturbation theory which 
requires small couplings. 

For the case of coupling between AF and SC order pa- 
rameters we can also construct a renormalization group. 
Within this approach we expect to find all the condi- 
tions for the occurrence of weak first order transitions. 
Detailed calculations and results will be published else- 
where. 



ONE-LOOP EFFECTIVE POTENTIALS AND 
RENORMALIZATION GROUP 



VI. CONCLUSIONS 



In the previous sections we have obtained weak quan- 
tum first order transitions from minima generated bal- 
ancing terms in the effective potential. This is possible 
when we compare terms of the same order and therefore 
the minima depend on the values of the couplings. Con- 
sidering, for example, the superconductor coupled to the 
electromagnetic field, the terms we have to balance are 
proportional to A and q 4 . New minima away from the 
origin are produced if 

A ~ q 4 (32) 

and this condition also leads to a relation for the Lon- 
don penetration depth and the coherence length A^ <C £. 
In this section we show how we can use renormalization 
group arguments to generalize the results for any small 
A and q and get rid of the condition (|32l . 

When deriving Eq. @ we have introduced a new pa- 
rameter M, associated with a renormalization mass and 
which was completely arbitrary. The results for the 
model must be the same for any chosen value of M. Now, 
if we prove that a variation in the value of M produces 
a correspondent variation in the coupling constant A, it's 
always possible to satisfy the condition (|3*2")l choosing a 
suitable value of M. To prove that, we have to consider 
Eq. J2J). We want to rewrite this equation with a value 
M', such that, 



^ = 2™ V+ 4^ +64^ 
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which is equivalent to Eq. (J2J with the reparametrization 



4! 



A 
4! 



3q 4 , M' 2 
In ■ 
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(33) 



We have studied the effects on a phase transition of 
coupling the order parameter of this phase to soft modes 
or to fluctuations of a competing phase. Since these ef- 
fects become more important at low temperatures we 
have mostly considered the zero temperature case. We 
have obtained the effects of quantum corrections due to 
this coupling on the quantum critical point associated 
with the relevant order parameter. These corrections 
were calculated using the effective potential method. 
They can produce radical modifications on the nature 
of the original phase transition. They can change the 
nature of the transition, from continuous to discontin- 
uous and also affect the ordered phase itself giving rise 
to an inhomogeneous phase with two values of the order 
parameter. 

Whenever the fluctuations change the nature of the 
transition to weak first order, we can study them using 
a scaling approach. We have obtained explicit results for 
one particular case but that should hold in general. As 
the system approaches the discontinuous transition from 
non-zero temperature, we can distinguish three different 
regimes. In the first, the high temperature regime, the 
scaling behavior is the same as approaching the quantum 
critical point of the uncoupled system, i.e., in the absence 
of additional fluctuations or soft modes. Further decreas- 
ing the temperature there is a non-universal regime which 
may depend on the particular dynamics of the original 
system and of the nature of the additional modes. Fi- 
nally at the lowest temperatures there is a regime which 
is characteristic of the first order nature of the zero tem- 
perature transition. At such low temperatures the corre- 
lation length saturates at a finite value. The thermody- 
namic properties have a contribution which is thermally 
activated. This can be traced to excitations which are 
gapped due to finite size effects as they occur in finite 
bubbles of the incipient phase. The correctness of this 
interpretation is clear from the fact that the gap for ex- 
citations is inversely related to the size of the bubbles 
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which in turn is of the order of the correlation length. 
Notice also that the former is related to the shift in the 
zero temperature phase transition with respect to that of 
the uncoupled system. 

The relevant exponents which determine the scaling 
behavior near the weak first order transition are those 
associated with the quantum critical point of the uncou- 
pled system, i.e., of the system in the absence of com- 
peting fluctuations or soft modes. In particular this is 
the case of the dynamic exponent. The scaling proper- 
ties of the new transition are related to the existence of 
an underlying second order instability which has become 



discontinuous due to the effects of soft modes or of fluc- 
tuations of a competing order parameter. 
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